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ABSTRACT 

This  report  documents  recent  efforts  in  digital  computer  modeling 
of  multipath  propagation  effects  on  a  modern  tracking  radar.  Specific 
questions  on  the  role  of  the  earth  dielectric  constant  and  surface 
roughness  are  considered.  Several  earth  models  are  also  considered, 
including  a  flat  smooth  surface,  smooth  4/3  radius  spherical  surface, 
and  a  4/3  radius  spherical  earth  with  rough  surface.  Two  basic 
scattering  models  are  considered,  specular  reflection  which  is  coherent 
in  nature  and  diffuse  reflection  which  is  noncoherent. 

This  study  uses  a  large  digital  phased  array  radar  simulation  which 
will  be  described  in  detail  in  a  future  report.  The  parts  of  signifi¬ 
cance  to  this  effort  are  described  for  completeness. 
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1.  Introduction 

The  main  purpose  of  this  report  is  to  document  recent  multi- 
path  propagation  studies,  particularly,  the  effects  of  this  phenomenon 
upon  the  angle  tracking  accuracy  of  a  modern  tactical  radar.  The  primary 
goal  is  to  provide  a  time  domain  digital  simulation  of  a  phased  array 
monopulse  tracking  radar.  To  properly  exercise  this  simulation,  all 
input  functions  must  be  carefully  modeled  to  provide  the  proper  temporal, 
amplitude,  and  phase  characteristics.  This  report  develops  the  multi- 
path  component  of  the  input  and  presents  results  obtained  during  verifi¬ 
cation  of  the  model. 

The  simulation  is  described  in  two  parts;  multipath  models,  and 
target  and  radar  models.  The  section  on  multipath  develops  the  several 
possible  components  of  the  multipath  signal  and  illustrates  the  effect 
of  each  component  upon  the  radar.  Finally,  the  composite  multipath  sig¬ 
nal  is  considered.  Sufficient  mathematical  expressions  are  ircluded  to 
aid  in  outlining  the  approach,  and  references  are  given  for  those 
interested  in  the  detailed  _(athematical  derivations.  Target  and  radar 
models  are  presented  only  to  show  how  the  tracking  errors  caused  by 
multipath  were  developed.  Details  of  the  radar  model  will  be  published 
in  a  later  report. 

This  report  addresses  multipath  errors  only  in  the  sin  p  coordinate, 
which  is  one  dimension  of  the  antenna  (sin  a,  sin  p)  steering  space  and 
corresponds  to  the  elevation  coordinate  only  when  sin  a  =  0.  For  other 
locations  multipath  errors  exist  in  the  sin  a  channel  also.  A  cursory 
examination  indicated  that  these  errors  would  be  quite  small.  Therefore, 
the  simulation  runs  were  normally  made  for  near  zero  sin  a.  A  very 
limited  number  of  runs  for  large  sin  a  steered  angle  produced  sin  p 
error  profiles  quite  similar  to  those  for  sin  a  =»  0  and  sin  a  multipath 

_3 

errors  that  were  considered  insignificant  (10  millisine  peaks).  There¬ 
fore,  multipath  error  in  the  sin  a  coordinate  was  not  addressed  in 
great  detail. 

The  results  (Paragraph  5)  indicate  the  sensitivity  of  multipath 
induced  tracking  error  to  two  classes  of  variables;  those  related  to  the 
input  and  those  related  to  the  radar.  Input  variables  of  interest  are 
the  relative  dielectric  constant  and  roughness  of  the  earth  surface,  the 
elevation  angle  of  the  target,  and  the  target  radial  velocity.  Variables 
associated  with  the  radar  are  antenna  patterns,  including  location  and 
magnitude  of  sidelobes  in  both  sum  and  difference,  antenna  height  above 
the  earth's  surface,  data  rate,  and  tracking  filter  bandwidth. 

Appendix  A  contains  the  details  of  the  computer  programs  used  for 
the  antenna  and  multipath  computations,  a  dictionary  of  variable  names, 
and  simple  program  outline.  These  programs  are  not  intended  as  final 
versions,  no  attempt  was  made  to  improve  their  efficiency. 
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2.  Multipath  Models 


The  basic  earth  and  scattering  models  will  be  developed  in 
this  section.  Only  five  of  the  many  combinations  will  be  considered: 


a)  Flat  smooth  earth,  specular  reflection 

b)  Spherical  smooth  earth,  specular  reflection 

c)  Spherical  rough  earth,  specular  reflection 

d)  Spherical  rough  earth,  diffuse  reflection 

e)  Spherical  rough  earth,  composite  reflections. 


The  flat  earth  model  will  be  considered  in  detail  because  of  its 
basic  nature.  Also,  its  solution  will  be  used  as  part  of  the  solution 
for  the  more  complex  cases. 


a ,  Flat  Smooth  Earth,  Specular  Reflection 

The  problem  of  signal  propagation  between  two  points  is 
illustrated  below.  Two  terminals  are  shown  t'c  different  heights  above 
a  flat,  smooth  ground  plane  characterized  by  a  relative  dielectric 
constant  (er) .  The  point  on  the  surface  where  the  reflection  occurs  is 

a  function  of  the  terminal  heights  and  the  distance  between.  The  problem 
is  reciprocal,  either  terminal  may  be  considered  the  transmission 
location. 


The  justification  for  this  rather  simple  geometrical  optics  model 
may  be  found  in  several  texts,  notably  Beckmann  and  Spizzichino^ .  To  a 


^3eckmann,  P. ,  and  Spizzichnio,  A.,  The  Scattering  of  Electromagnetic 
Waves  from  Rough  Surfaces,  MacMillan  Co.,  New  York,  1963. 
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first  approximation,  the  reflection  point  may  be  considered  as  the 
center  of  the  first  Fresnel  zone,  i.e.,  the  set  of  points  on  the  surface 
with  path  length  difference  less  than  a  half  wavelength.  Geometrical 
optics  will  be  used  throughout  this  report  to  determine  the  location 
of  the  reflection  point. 

The  amplitude  of  the  reflected  energy  is  determined  by  the 
reflection  coefficient  at  the  earth's  surface.  For  the  case  of  vertical 
polarization  of  the  transmitted  and  received  signals,  the  reflection 

2 

coefficient  given  by  Kerr  is 


where 

K^  =  electromagnetic  propagation  factor  (complex) 

Kg  =  freespace  propagation  factor  (complex) 

f2  =  grazing  angle 

Tv  =  complex  reflection  coefficient 

Pv  =  magnitude  of  reflection  coefficient 

and 

K  =  a  +  jwe 

a  =  conductivity  of  medium 
a)  =  radian  frequency  of  energy 
e  =  permittivity  of  medium. 
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Kerr,  D.  E.,  Propagation  of  Short  Radio  Waves,  McGraw-Hill  Book  Co., 
New  York,  1951. 
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For  most  soil  conditions  under  consideration  e  >  >o  so  that  K  =  jwc 
and  the  above  expression  becomes 


°v  = 


er  sin  ^ 


:ZH 


cos  \j»2 


&r  sin  ^2  +  /  tr  “  cos 


Therefore , 
and  has  no 
•^2  becomes 

polarizing 
Above  this 


the  reflection  coefficient  is  no  longer  a  complex  quantity 
effect  on  the  phase  of  the  reflected  signal.  However,  when 
less  than  a  critical  angle,  (the  complement  of  Brewster's 

angle),  the  reflected  signal  has  a  phase  shift  of  it  radians, 
angle  no  phase  shift  occurs. 


To  illustrate  how  the  signals  must  be  combined  a  vector  diagram  will 
be  useful.  Assume  that  the  direct  path  energy  is  taken  as  the  reference. 
Its  amplitude  and  phase  depend ,  in  general,  upon  such  things  as  transmitted 


power,  receiver  antenna  area,  losses,  and  propagation  path  length  and 
characteristics.  For  this  discussion  it  may  be  taken  as  unity  amplitude 
and  zero  phase  angle.  The  indirect  signal  is  scaled  by  the  reflection 
coefficient  magnitude  and  has  a  phase  a  determined  by  the  path  length 
difference  and  the  possible  n  reflection  phase.  Then 

“  =  f  (RID  -  Rd)  £or  *2  >  *c 

and 

a  -  f  (rid  -  rd)  +  *  £or  H  -  *C 


wbeie  \  =  wavelength. 


A 


This  method  of  calculating  a  requires  determination  of  and 

quite  accurately.  This  is  well  within  the  capability  of  the  CDC  6600 
computer  used  for  this  study,  but  problems  have  been  encountered  with 
other  equipment.  One  method  used  to  overcome  this  problem  is  to  assume 
^2  53  ^i  >  thus  allowing  a  simple  calculation: 

a  =  f  (2hl  sin  *l) 

where  h^  is  the  height  of  terminal  1,  and  \|^  is  the  elevation  angle  of 
terminal  2  from  terminal  1  (see  sketch  page  2). 

b.  Spherical  Smooth  Earth,  Specular  Reflection 

The  geometry  of  the  spherical  earth  multipath  problem 
with  exaggerated  dimensions  is  shown  in  the  following  sketch.  In  this 
model  the  geometrical  optics  reflection  is  with  respect  to  the  earth 
tangent  at  the  point  of  reflection,  so  that  in  general  the  multipath 
points  for  a  flat  earth  and  a  curved  earth  model  are  not  the  same. 
Several  methods  have  been  proposed  for  the  calculation  of  the  reflection 

location,  and  the  best  for  use  in  this  case  was  found  in  Kerr^.  Once 
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Kerr,  op.  cit. 
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fha  geometrical  problem  is  solved  the  reflection  coefficient  for  vertical 
polarization  may  be  calculated  from 


sv 


=  0  D 
v 


where  pv  is  the  flat  earth  vertical  polarization  reflection  coefficient, 
and  D  is  the  divergence  factor  caused  by  spherical  reflection 


D  = 


r 


All  other  parts  of  the  calculation  such  as  a^^nna  scaling  and  phase 
calculations  are  handled  exactly  ljce  the  flJ*  earth  cas*’.  Derivations 
of  the  divergence  factor  D  and  the  tv*  i-  of  the  refl  rtion  location 

are  given  in  Kerr^ . 


c .  Spher leal  Rough  Ea;  _  Specular  Reflection 

When  , v.irtce  roughnefs  .  s  added  to  the  model  several  new 
parameters  become  involved,  and  the  character  of  che  model  changes.  The 
mod^.s  described  so  Jar  have  fc..'n  deterministic  in  nature,  but  statisti¬ 
cal  methods  must  be  used  _o  defer ‘be  zb*.  distribution  of  heights  and 
spacings  of  the  su.jb?'\  ,irotr«.si‘  ..  and  •  cattering  from  this 
surface . 


In  the  case  of  a  smooth  surface  al^  the  indirect  field  is  scattered 
specularly  by  the  Fresnel  zone  described  previously.  As  the  surface  is 
broken  up  by  roughness  this  field  decomposes  into  two  components .  One 
is  t.ae  specular  component,  existing  in  the  specular  direction  as  defined 
by  the  smooth  s  to  face  case,  If«  charxr ter? '-tics  will  be  described  in 
this  paragraph.  The  other  component  is  diffuse  or  noncoherei t  and 
exists  in  all  directions;  it  w^’l  be  described  in  the  next  paragraph. 

If  the  surface  is  allowed  to  have  varying  vertical  heights  about 
the  smooth  surnce  as  a  mean,  then  the  waves  reflected  from  two 
adjacent  poj  ;  within  the  first  Fresnel  zone  of  height  difference  Ah 
«ill  differ  ?n  phase  by 


4flAh  . 
a  sir.  '^2 


.err,  loo,  cit. 
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Further,  if  it  is  assumed  that  the  surface  heights  are  normally 
distributed  with  standard  deviation  Ah  and  zero  mean,  then  according 

to  Beckmann^  the  distribution  of  the  scattered  energy  consists  of  two 
components,  a  fixed  component  and  a  fluctuating  component  that  is  Hoyt 
distributed  (orthogonal  normal  variates  with  zero  means  and  unequal 
variances) .  The  composite  scattered  field  has  a  mean  square  value  of 


1 

■? 

i 


2 


0 


s 


e-(A4>) 
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This  is  then  used  with  the  reflection  coefficient  of  (b)  to  give  the 
rough  earth  specular  reflection  coefficient 


For  short  tranrmission  paths  over  land  the  coefficient  may  be  taken  as 
a  constant  given  by  the  mean  value  expression.  This  was  done  in  the 
study.  Also,  the  adjacent  surface  features  were  assumed  independent, 
and  calculated  and  used  as  indicated.  However,  other  very  important 

parameters  are  the  correlation  distance  and  the  surface  autocorrelation 
function.  This  surface  correlation  might  be  introduced  in  a  manner 

similar  to  that  used  by  Pope  and  Wood^  to  introduce  correlation  into 
simulated  ground  clutter  data. 

d .  Spherical  Rough  Earth,  Diffuse  Reflection 

The  diffuse  component  of  the  rough  earth  scattering 
emanates  from  an  area  of  the  surface  which  is  in  general  much  larger 

than  the  first  Fresnal  zone.  Beckmann  defines  this  area  in  terms  of 
the  terminal  geometry  and  mean  surface  slope.  Once  this  area  is  defined, 
calculation  of  the  diffuse  scattering  coe  fficient,  oD,  involves  integration 
over  the  surface,  the  radar  antenna  pattern,  and  the  target  scattering 
characteristics.  In  order  to  simplify  the  model  for  this  study,  the 


Beckmann,  loc.  cit. 

^Pope,  B.  E.,  and  Wood,  W.  E.,  Digital  Simulation  of  a  Phase  Monopulse 
Tracking  Radar,  U.S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama, 
Report  No.  RE-TM-70-2,  28  January  1970  (Unclassified). 

^Beckmann,  loc.  cit. 
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RMS  diffuse  scattering  coefficient  was  defined  as  in  the  following 
diagram.  This  may  be  compared  with  the  generalized  results  given  by 

O 

Beckmann  on  page  340.  The  distribution  was  assumed  Rayleigh,  after 

Q 

Beckmann  . 

The  point  of  origin  of  this  scattering  must  also  be  determined.  In 
most  cases  of  interest  in  this  study  the  diffuse  return  may  be  assumed 
to  originate  at  the  specular  point.  This  may  be  invalid  for  elevation 
angles  less  than  a  few  tenths  of  a  degree.  However,  simulation  runs 
with  the  diffuse  scattering  from  other  locations  produced  similar  results. 
Therefore  the  specular  and  diffuse  scattering  components  were  assumed 
co-located  in  this  study. 


In  this  case  the  scattering  coefficient  becomes 


prsv  "  psvpDD 

The  divergence  factor  D  may  be  removed  from  this  expression  according 

to  Beckmann^,  because  it  has  little  effect  on  the  diffuse  scattered 
signal . 
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Beckmann, 

Beckmann, 

Beckmann, 


loc.  cit. 
Ibid . 
loc.  cit. 


8 


'  *"*''*■  VI"  ?.TV  ” 


e .  Spherical  Rough  Earth,  Combined  Reflections 

Rough  surfaces  may  be  taken  as  those  which  satisfy  the 
Rayleigh  criterion, 


A4>  <  ^ 


This  may  also  be  expressed  in  terms  of  the  grazing  angle  and  RMS  surface 
height  as 


sl”  +2  <8Zh 


Thus,  it  can  be  seen  that  a  surface  may  be  rough  in  the  Rayleigh  sense 
in  one  grazing  angle  region  and  smooth  in  another.  This  could  be  used 
to  define  the  angular  regions  of  importance  of  various  types  of 
scattering. 


The  diffuse  and  specular  models  may  be  used  together  in  a  composite 
model.  The  same  aquations  are  used  and  the  signals  combined  as  indicated 
in  the  vector  diagram  below. 
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3.  Radar  and  Target  Models 
a .  Radar  Model 


The  radar  model  used  for  this  study  is  described  below 
in  block  diagram  form. 


The  basic  antenna  model  (Antenn  is  a  rectangular  array  space  fed  by 
four  horns.  On  receive,  these  ti,  ..a.  signals  are  combined  to  form  ampli¬ 
tude  monopulse  sum  and  two  difference  voltages.  These  three  signals  are 
split  into  in-phase  and  quadrature  components  by  the  down  conversion/mixer 
which  is  assumed  ideal,  as  was  the  monopulse  network.  The  phase  of  this 
local  oscillator  is  selected  randomly  over  two  X  radians.  After  analog- 
to-digital  (A/D)  conversion  an  angle  measurement  is  made  by  the  following 
algorithm: 

(SI  X  Dl)+  (S0 

Angle  Measurement  =  ~ 

(SI  *  SI  +  SQ 


xDo) 


X  s, 


q) 


X  Slope 


where 


Sj.  =  Sum  in-phase  component 

Sq  =  Sum  quadrature  component 

Dj.  =  Difference  in  phase  component 

Dq  =  Difference  quadrature  component 

Slope  =  Monopulse  slope  factor  to  convert  error 
voltage  to  angle. 

This  measured  angle  is  then  smoothed  by  a  critically  damped  GH  filter 
and  a  prediction  made  for  the  next  beam  position.  The  filter  beta  was 
maintained  at  0.8  and  a  sampling  race  of  10  samples  per  second  was  used 
unless  otherwise  stated. 
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b.  Target  Model 

For  simplicity  a  constant  cross  section  target  of  1 
square  meter  was  flown  radially  outward  from  the  radar. 


Z 


-Y 


Initial  values: 

XQ  =  200  m 
Yq  =  200  m 
Zq  =  50  m 
Xq  =  141  m/ sec 
Yq  =  141  m/sec 


4.  Stability  Requirements 

During  the  investigation  it  was  noted  that  some  combinations 
of  data  rate,  target  velocity,  and  track  filter  bandwidth  produced 
unstable  angle  tracking  and  in  some  cases  complete  loss  of  track. 
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Figure  1  shows  a  portion  of  the  error  profile  where  unstable  operation 
was  encountered  but  stable  track  recovered.  The  data  rate-velocity- 
bandwidth  space  was  mapped  to  locate  the  boundary  between  the  stable 
and  unstable  regions  for  specular  reflections  only.  The  instability  is 
observed  as  an  oscillation  as  in  Figure  1,  and  loss  of  track  arbitrarily 
defined  as  an  error  in  excess  of  30  millisines.  The  space  thus  defined 
is  shown  in  Figure  2,  with  the  region  of  stable  operation  located  below 
the  pictured  surface  (high  velocity  and  low  sample  rates  are  highly 
stable,  low  velocity  and  high  data  rate  tend  toward  instability). 

To  explain  this  behavior,  a  stability  analysis  was  undertaken. 
Since  the  combination  of  direct  and  indirect  signals  produced  a  time 
varying  nonlinearity  the  analysis  was  done  first  for  a  no-multipath 
environment.  This  analysis  is  presented  in  Appendix  B  and  outlined 
below. 

FORWARD 

TRANSFER  FUNCTION 


INPUT 


% 

Q 

OUTPUT 

( 

1 _ 

H 

FEEDBACK 

The  forward  transfer  function  may  be  considered  as  two  separate  blocks 
for  clarity. 


The  previous  sketch  represents  actual  system  parameters  where: 

^  =  Difference/ sum  voltage  ratio 

AP  «  Input  angle  error 

6  ■  Measured  angle  error 
m 
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|  =  kt  A? 

P  =  k0 
Hm  2  l  £ 


For  Antenna  1  and  matched  conditions 


k^  ~  82.5  volt/volt/radian 
~  0.01208  radian/ volt/volt 


For  the  actual  analysis  of  the  system  these  two  may  be  combined  into  one 
transfer  block  with  gain  K  given  by  k^k2*  Note  that  for  normal  or 

"matched"  operation  K  =  1.  This  transfer  function  is  shown  below: 


The  analysis  using  the  forward  transfer  function  produced  the  following 
stability  criteria  as  a  funtion  of  beta,  the  bandwidth  factor  for  the 
feedback  path  (see  the  following  sketch). 

The  analysis  proved  that  as  long  as  the  system  transfer  function 
is  linear  and  matched  values  are  used  (k^  =  l/k2)  the  system  is  stable. 

However,  the  processor  gain  may  be  improperly  chosen  so  that  instability 
occurs . 

A  few  comments  on  A/Z  curves  are  appropriate.  As  shown  in  Appendix 
B,  these  curves  are  nonlinear  and  in  some  cases  may  be  linearized  over 
only  a  limited  region.  Selection  of  a  linear  function  for  the  inverse, 
k2,  is  always  somewhat  in  error.  Therefore  the  idea  of  matched  antenna 

and  processor  is  only  an  ideal. 

When  multipath  propagation  is  added  to  the  problem,  the  forward 
transfer  function  becomes  time-varying  nonlinear.  However,  the  remainder 
of  the  system  is  unchanged,  so  that  the  stability  criteria  derived  for 
the  no-multipath  case  still  apply.  Referring  to  Figure  11,  the  instabil-  • 
ity  only  occurs  at  the  multipath  peaks  at  0.8  and  0.4  degree  elevation 
angles.  The  limitation  to  the  two  locations  may  be  explained  by 
reference  delta/sum  (D/S)  curves  produced  at  these  points. 
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D/S  curves  and  transfer  functions  were  generated  for  the  error 
peak  located  at  0.8-degree  elevation  to  show  how  the  slope  stability 
criteria  may  be  applied.  The  D/S  curves  are  given  in  Figures  3,  4,  5, 
and  6.  The  first  and  last  of  these  figures  are  for  stable  locations, 
while  the  other  two  are  on  the  unstable  peak.  Each  curve  was  generated 
by  placing  a  stationary  target  at  the  elevation  angle  indicated,  and 
the  antenna  swept  to  produce  the  D/S  curve.  These  figures  may  be 
compared  with  Figur«.  11  to  determine  their  relationship  with  error 
direction  and  magnitude. 

Transfer  functions  of  the  form  Pm  =  KA(3  were  then  developed  from 

the  D/S  curves  and  k2  =  0.01208  radian/volt/volt.  These  are  shown  in 

Figures  7,  8,  9,  and  10  which  correspond  to  the  preceding  D/S  curves. 
Stability  limits  for  selected  values  of  beta  from  the  stability 
analysis  are  shown  in  Figure  8. 

The  difference  between  complete  loss  of  track  as  defined  in  this 
case  and  unstable  track  with  recovery  is  a  function  of  the  data  rate 
and  target  velocity.  Stated  qualitatively  this  reduces  to  the  following: 
the  more  samples  taken  in  the  unstable  region,  the  greater  the  effect 
on  angle  tracking.  Therefore,  low  velocities  and  high  data  rates  allow 
many  measurements,  while  high  velocities  and  lower  data  rates  reduce 
the  number  of  measurements  taken  while  the  system  is  prone  to  instability. 
(This  agrees  with  the  volume  in  Figure  2.)  The  amount  that  K  exceeds 
the  stability  requirement  also  effects  the  number  of  samples  needed  for 
loss  of  track. 
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5.  Results 

Computer  runs  were  organized  to  determine  the  sensitivity  of 
multipath  error  to  target  altitude,  target  velocity,  earth  dielectric 
constant  and  roughness,  antenna  height  and  pattern,  and  the  path  length 
difference  algorithm.  For  convenience  the  target  was  started  close  to 
the  radar  and  flown  radially  outward.  A  few  selected  cases  were  run  for 
incoming  targets  and  the  results  were  almost  identical. 

It  was  determined  early  in  the  effort  that  the  only  two  geometrical 
variables  of  interest  were  the  radar  antenna  height  and  the  elevation 
angle  of  the  target.  Thu3  target  altitude  and  range  do  not  appear  on 
the  plots . 

Figure  11  shows  the  low  sensitivity  of  multipath  errors  to  earth 
relative  dielectric  constant. 

Figure  12  shows  a  similar  but  more  pronounced  effect  when  surface 
roughness  is  considered  in  a  specular  model.  Larger  values  of  RMS 
surface  height  tend  to  reduce  the  specular  component,  but  increase  the 
diffuse  component  as  shown  in  Figure  13.  Note  that  Figure  13  illustrates 
the  envelope  of  the  error,  which  appeared  as  a  noise  voltage.  This  was 
expected  because  the  phase  angle  was  drawn  from  a  random  uniform 
population.  When  both  specular  and  diffuse  scattering  are  present  the 
result  appears  as  a  superposition  of  the  incoherent  diffuse  component 
and  the  highly  coherent  specular  component.  The  envelope  of  the  errors 
are  shown  in  Figure  14.  It  can  be  seen  that  the  random  diffuse  errors 
are  distributed  about  the  specular  error  as  a  mean  value. 

A  much  larger  variation  in  results  occurred  when  changes  in  antenna 
height  and  pattern  were  allowed.  Figures  15  and  16  illustrate  the 
antennas  used,  the  major  difference  being  in  the  width  of  difference 
pattern  major  lobes  and  the  location  of  nulls.  Antenna  1  has  coincident 
sum  and  difference  nulls  so  that  narrow  difference  lobes  result.  Antenna 
2  has  broad  difference  lobes.  Both  sum  beams  have  a  3  dB  beamwidth  of 
1.6  degrees  in  azimuth  and  elevation.  The  two  antenna  heights  were 
chosen  for  comparison  with  other  published  results.  Figures  17  and  18 
show  the  type  results  encountered.  Changing  the  antenna  height  relocates 
the  maximum  multipath  errors  in  the  region  below  2  degrees.  Not  only 
are  the  locations  changed,  but  the  maximum  values  change.  Changing  the 
antenna  model  (Figure  18)  also  produced  drastic  changes.  Note  that  this 
run  is  plotted  on  a  different  scale,  because  the  errors  were  considerably 
larger  than  previous  values.  Several  variants  of  both  antennas  were  run 
at  both  heights  with  similar  results:  the  location  and  magnitude  of  the 
peaks  were  in  general  unpredictable,  but  no  errors  above  18  to  20 
millisines  were  observed. 

The  approximate  algorithm  for  calculating  path  length  difference 
mentioned  in  Paragraph  4  was  also  run  for  comparison.  Figures  19  and  20 
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show  the  results.  In  the  first  case  Antenna  1  was  used  at  a  4  meter 
altitude,  in  the  second  case  Antenna  2  was  used  at  a  3  meter  altitude. 

In  general  this  algorithm  produced  results  similar  to  the  exact  solution 
with  slightly  different  peak  magnitudes  and  locations. 

In  addition,  some  runs  were  made  with  the  antenna  constrained  to 
operate  only  above  certain  elevation  angle  limits,  sometimes  called 
super-elevation  tracking.  The  multipath  errors  were  reduced  in  magnitude 
but  not  eliminated  entirely. 


? 
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6.  Conclusions 


This  study  has  shown  that  multipath  errors  can  be  significant 
in  a  narrow  beam  tracking  radar.  These  errors  are  almost  independent 
of  the  relative  earth  dielectric  constant.  Surface  roughness  produces 
a  more  pronounced  effect  on  both  specular  and  diffuse  returns.  The 
location  and  magnitude  of  the  major  multipath  errors  are  extremely 
sensitive  to  variations  in  antenna  height  and  difference  pattern  structure, 
especially  in  the  region  of  the  peak  lobe  and  first  null.  With  the 
tracking  filter  in  the  loop  most  of  the  errors  were  below  10  millisines, 
and  all  were  below  16  millisines.  Some  combinations  of  target  velocity, 
angle  tracking  filter  characteristics,  and  data  rate  produced  unstable 
track  situations. 
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Figure  5.  A/Z  Curve  with  Multipath,  Elevation 
Angle  s  0.75  Degree,  Antenna  1 
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Figure  6.  A/Z  Curve  with  Multipath,  Elevation 
Angle  =  0.70  Degree,  Antenna  1 
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Figure  9.  Radar  Transfer  Function,  Elevation  =  0.75  Degree, 
Antenna  1 
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Appendix  A 

PROGRAM  LISTING  AND  DEFINITION  OF  VARIABLES 


This  appendix  contains  a  listing  of  the  multipath  and  antenna 
subroutines  along  with  a  variable  list  and  comments. 

The  operation  may  be  divided  by  line  numbers  as  follows: 

SUBROUTINE  MULT 

100  -  160  Dimension,  complex  and  common  statements 

180  -  320  Determination  of  multipath  point  over  spherical 

earth  and  conversion  to  radar  coordinates 

340  -  400  Determine  phase  relationship  of  direct  and  indirect 

signals 

420  -  610  Determine  reflection  coefficient  for  specular  and 

diffuse  returns. 

630  Calculate  antenna  response  to  direct  signal 

650  Calculate  antenna  response  to  indirect  signals 

770  -  820  Combine  direct  and  indirect  signals. 

Dictionary  of  selected  variables: 

AE  Radius  of  4/3  earth,  8.5  E  +  06  meters 

ALPM  Phase  difference  between  direct  and  indirect  signals 

D  Spherical  earth  divergence  factor 

DFA  Complex  antenna  response  in  azimuth  difference 

channel 

DFE  Complex  antenna  response  in  elevation  difference 

channel 

DXR 

DYR  Target  velocities  in  array  face  rectangular  coordinate 

DZR 

GAIN  Peak  sum  voltage  gain 

HA  Antenna  height,  meters 

HT  Target  Height,  meters 

PI  3.14159 

PSIC  Critical  reflection  angle,  radians 

PS12  Grazing  angle  at  multipath  point,  radians 


Preceding  page  blank 
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REFCD 

REFCS 


SIGMH 

SMUET 

SMULP 


SXALP 

SXBET 

WAVLN 


Target  range,  meters 

Flat  earth  vertical  polarization  reflection 
coefficient 

Diffuse  reflection  scattering  factor 
Specular  reflection  scattering  factor 

True  target  angles  in  radar  sine  space 

RMS  surface  roughness  height 

Multipath  point  angles  in  radar  sine  space 

Complex  antenna  response  in  sum  channel 

Predicted  target  position  in  radar  sine  space 

Radar  wavelength 


YM  Multipath  point  in  rectangular  coordinates 

ZM 

XN  Number  of  elements  per  row  or  column  of  array 

ZN1  Real  part  complex  relative  earth  dielectric  constant 

The  listing  of  Subroutine  MULT  gives  the  exact  method  of  calculating 
path  length  difference.  To  compare  this  with  the  approximate  solution, 
lines  2770  through  2980  were  changed  or  deleted  to  reflect  the 
equations  given  at  the  end  of  Paragraph  2. a. 

Subroutine  ANTENA. 


930  -  960 


initial  conditions 


970  -  1100 


Calculate  antenna  response 


Dictionary  of  selected  variables: 


Three  dB  sum  becmwidth 
Azimuth  difference  channel  voltage 
Elevation  difference  channel  voltage 
Normalization  factor 


Responses  at  each  of  the  four  illumination  horns 


Measurement  angles  in  sine  space 

Beam  pointing  angles  in  sine  space 

Sum  channel  voltage  response 
Wavelength  of  energy 

Number  of  elements  in  each  row  and  column  of  the 
rectangular  array 

Two  antenna  models  are  referenced  in  this  report.  The  listed 
version  of  Subroutine  ANTENA  calculates  the  response  of  Antenna  Modal  1. 
Antenna  Model  2  was  implemented  by  substitution  of  the  following  for 
line  1080: 


SALP 

SBET 

SNALP 

SNBET 

SUM 

WAVLN 

XN 


XI  m  SBET  -  SNBET 
X2  *=  ABS  (XI) 

DX  =  X2  *  180. PI 

IF  (DX.GE.0..AND.DX.LE.3.)  DE  =  0.707  *  SIN(180./3.*X1) 

IF  (DX.GT.3..AND.DX.LE.4.)  DE  =  0.063  *  SIN  [-180.  *  (Xl-3.  *  PI/180.)] 
IF  (DX.GT.4.)  DE  =  0.01  *  SIN  [180.  *  (Xl-4.  *  PI/180.)] 

DFE  =  CMPLX  (DE,0.) 


SUBROUTINE  MULT (SNALP.SNBET .SALP.SBET .R.SUM.DFA.DFE.NN) .RETURNS 
COMPLEX  XIDD.XIOS.SUM.DFA.DFE.SUMD.DFAO.DFED 
COMPLEX  SUMIO.OFAIO.OFEID 

COMMON/MA1N/XP.YR.ZR.XT.  YT  .ZT  .WAVLN.T  .HA.DXR.DYP.DZR.XL 
NAMELIST/NAMMU/ZN1.SIGMH.XSPEC.XDIFF.GMA.8ETAO 
IF(R.LT.l.)GO  TO  100 

DATA  PI. C.AE/3. 1415926536. 2. 997925E*08. 8. 5E*06/ 

OETERMINE  LOCATION  OF  MULTIPATH  POINT  IN  TARGET  CO-ORDINATFS 

Z=ZT*AE»HA  SHT=S0RT (XL**2+Z**2> -AE 

S=AE»ACOS (Z/SQRT (XL**2*Z**2) ) 

P=2.#S0RT (AE/3.#(HT*HA) ♦ ( (S/2. )«*2)/3. ) 
PHI=AC0S(2.*AE*S*(HA-HT)/P**3> 

Rl-S/2.*P*C0S((PHI*PI)/3.)  SR2=S-R1 
ZL=AE*CCS(R1/AE)  SXLM=SQRT (AE**2-ZL**2) 

ZMsZL-AE-HA 

XM=XLM*C0S(ATAN2(YT.XT))  SYM=XLM*SIN(ATAN2(YT»XT> ) 

IF((AE**2)/ZL.GE.(AE*HA)>  RETURN  M 

CONVERT  TO  RADAR  CO-ORDINATES 

CALL  CNVRT (XM.YM.ZMtA.B.D.SMULP.SMUET.GMA) 

DETERMINE  PATH  LENGTH  DIFFERENCE  AND  CRITICAL  ANGLE 

RID=SQRT (XM**2+YM»*2»ZM**2> ♦SORT { (XT-XM) **2* ( YT-YM) **2* (ZT-ZM) * 
ALPM=2.*PI* (RID-R) /WAVLN 

PSI2=ATAN(XM/ZL) ♦AC0S( (ABS (XM-XT) l/SORT ( (XM-XT>**2MZM-ZT)**2) ) 
PSIC*  ASIN  (SORT  (l./(ZNl+l.))) 

*F (PSI2.LE.PSIC) ALPMrALPM.Pl 

CALCULATE  REFLECTION  COEFF  AT  MULTIPATH  POINT 

TERM1=ZN1*SIN(PSI2) 

TERM2=SQRT (ZN1- (COS (PSI2) **2) ) 

REFC= ABS ( ( TERM 1 -TERM2  J / ( TERM 1+TERM2 ) ) 

0=1. /SORT (1 •♦4.*R1*R2/ (AE*S*SIN(2»*PSI2) ) ) 
DELPHI=4.*SIGMH*SIN(PSI2)/WAVLN 
REFCS=SQRT(EXP(-(DELPHI*PI)**2) ) 

PICK  RAYLEIGH  DISTRIBUTED  SAMPLE 

SIGMD=0.35/0.25*DELPHI 

IF (DELPHI .GT.O .25) SIGMD=0 .35 

SIGMD=S1GMD/SQRT(2.) 

CALL  NORMAL (SIGMD.O.. VI) 

CALL  NORMAL (SIGMD.O. »V2) 

REFCD=S0RT(V1**2+V2**2) 

ALPM0=RANF(U)*2.*PI 

XIDS=REFC*D*REFCS#CMPLX (COS (ALPM) .SIN(ALPM) ) 
XIOD=REFCD*REFC*CMPLX(COS(ALPMD) .SIN(ALPMD) ) 

:  CALCULATE  OIRECT  PATH  ANTENNA  GAIN 


(M)MULOOlOO 
MUL001 10 
MUL00120 
MULOO) 30 
MUL00140 
MUL00150 
MUL00160 
MUL00170 
MUL00180 
MUL00190 
MUL00200 
MUL00210 
MUL00220 
MUL00230 
MUL00240 
MUL00250 
MUL00260 
MUL00270 
MUL00280 
MUL00290 
MUL00300 
MUL00310 
MUL00320 ’ 
MUL00330 
MUL00340 
MUL00350 
•21MUL00360 
MUL00370 
MUL00380 
MUL00390 
MUL00400 
MUL00410 
MUL00420 
MUL00430 
MUL00440 
MUL00450 
MUL00460 
MUL00470 
MUL00480 
MUL00490 
MUL00500 
MUL00510 
MUL00520 
MUL00530 
MUL00S40 
MUL00550 
MUL00560 
MUL00570 
MUL00530 
MUL00590 
MUL00600 
MUL00610 
MUL00620 
MUL00630 
MUL00640 
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CALL  ANTFN4 ( SNALP « SNBf  T , SALP.SBFT , SUM, OF A .DFE » WAVLN.NN) 

SUMO=SUM 

OF  AO=OF A 

0FE0=0FE 

CALCULATE  INDIRECT  PATH  ANTENNA  GAIN 

CALL  ANTENA (SMULP.SMUET .SALP.SBET  .SUM.DFA.DFE.WAVLN.NN) 

SUMID=SUM 

OFEID=DFE 

0FAID=DFA 

COMBINE  DIRECT  AND  INDIRECT  SIGNALS 

5UM=SUMD*XSP£C»SUMID*XIDS+XDIFF*SUMID*XIDD 

TFA=DFAO*XSPEC*DFAID*XIDS+XDIFF*DFAID#XIDD 

OFE=OFED«XSPEC*OFEID«XIDS+XDIFF»DFEID*XIDD 

RETURN 

100  RE AD (5.NAMMU) 

VIRITF  (6.NAMMU) 

RETURN 

END 


MUL006S0 

MUL00660 

MUL00670 

MUL006B0 

MIJL00690 

MUL00700 

MUL00710 

MUL00720 

MUL00730 

MUL00740 

MUL007S0 

MUL00760 

MIJL00770 

MUL00780 

MUL00790 

MUL00800 

MUL00810 

MUL00820 

MUL00830 

MUL00840 

MUL00850 

MUL00860 


SUBROUTINE  ANTENA (SNALP. SNBET.SALP.SRET. SUM.DFA.DFE.WAVLN.NN) 

This  SUBROUTINE  DETERMINES  the  response  in  free  SPACE  OF  SINX/X  PAT 
SNALP  AND  SNBET  ARE  THF  BORESIGHT  ANGLES  IN  SINES 
SALP  AND  SHET  ARE  THE  MEASUREMENT  POINTS 

COMPLEX  UAHE,  UHHE.UAHA.U8HA.SUMyDFA.DFE 
NAMELIST/NAMAN/UAHA.UAHE.UBHE.UBHA.S 
IF(NN.EQ.O)GO  TO  100 

DATA  PI, C.AE/3. 1415926536. 2. 997925E*08.8.5E*06/ 

XN=NN 

P1=PI*S»XN*1./WAVLN 

P2=P1/XN 

P3=*AVLN/(?.*XN*S> 

FO=I./XN 

F1=SIN(P1*(SNALP-SALP*P3) )/SIN(P2*(SNALP-SALP+P3) )*FO 
F2=SIN(P1*(SNALP-SALP-P3) )/SIN(P2*(SNALP-SALP-P3) )*FO 
F3=SIN(P1*(SN8ET-SRET*P3))/STN(P2*(SNRET-SBET*P3) )*FO 
F4  =  SIN (PI  * (SNBET-SBET-P3) ) /SIN (P2« (SNBET-SBET-P3) ) *FO 
SUM=CMPLX(. 5* (Fl*F2)*(F3+F4)/0. 81069946,0.) 

OF  A" (F?-UBHA*F l)*(F3*F4)/( .81069946) *0.5 

OFE= (F4-UHHE*F3) *  IF  1 *F 2)  /  (0.81069946) *.5 

OFA=UAHA*DFA 

0FE=UAHF*0FE 

RETURN 

100  READ (S.NAMAN) 

WRITE (6.NAMAN) 

PE TURN 
END 


MIJL00870 
MUL00880 
MUL00890 
MUL00900 
MUL00910 
MUL00920 
MUL00930 
MUL00940 
MUL00950 
MUL00960 
MUL00970 
MUL00980 
MUL00990 
MUL01000 
MUL0I010 
MUL01020 
MUL01030 
MUL01040 
MUL01050 
MUL01060 
MUL01070 
MUL01080 
MUL01090 
MUL01100 
MULOlllO 
MUL01 120 
MUL01 130 
MUL01 140 
MUL01150 
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Appendix  B 

LINEAR  STABILITY  ANALYSIS 


The  angle  tracking  loop  of  this  radar  model  may  be  analyzed  as 
below . 


h  ^t-0 


te  + 


, - H  | - 

r - *1 

“2 

- 1 

A 

0 

A 

^ — 1 

H 

^ _ 1 

where 

G  =  forward  transfer  function  describing  the  radar  angle 
measurement  equations 

H  =  feedback  loop  gain,  here  a  critically  damped  gh  filter, 
ftp  =  true  target  position  (SNBET) 

A 

P  =  predicted  target  position  for  next  measurement  (SHBET) 

=  measured  angle  (SMBET) 

P^  =  measured  monopulse  error  signal  (VM) . 

The  names  enclosed  in  parenthesis  refer  to  varable  names  in  the 
computer  program  listed  in  Appendix  A .  Following  the  computer  listing, 
the  forward  block  may  be  modeled  as  p^  =  K  (p,j,  -  p),  or  in  sampled  data 

transfer  function  notation  G^(z)  s  K.  Thus  the  forward  loop  is  consid¬ 
ered  as  an  impulse  of  weight  K  occurring  at  time  n.  Similarly  G2(z)  =  1. 

The  transfer  function  of  the  feedback  loop  may  be  derived  from  the  three 
basic  gh  filter  equations: 


X  =  X  +  g  (x  -  x  ) 
n  n  °  V  n  n/ 

X  =  X  +  £  (X  -  £  ) 
n  n-1  T  V  n  n  / 


X  =  X  +  X  T 
rH-1  n 


Preceding  page  blank 


43 


where 


Xn  =  smooth  position  at  time  n 

A 

=  predicted  position  at  time  n 

Xn  =  raw  measured  position 

T  =  Sample  interval 
Xn  =  smoothed  velocity  at  time  n 


g  and  h  =  filter  parameters. 


These  three  equations  reduce  to: 


•Vi  -  Xn  <8  +  h)  -  g  X^  +  (2-g-h)  -  <1  -  g)  X^ 


If  the  output  is  considered  X  .  and  the  input  is  X  then  the  z 

nr  i  n 

transform  is 

H(z)  =  -S2 - 

z  -  (2-g-h)z  +  (1  -  g) 


Because  the  predicted  value  is  always  used  for  the  next  measurement 
calculation  a  one  period  delay  must  be  included.  Therefore  H(z)  becomes 


H(z)  =  — -  S - 

-  (2-g-h)z  +  (1  -  g) 


This  is  the  expression  used  in  the  analysis. 

To  simplify  the  analysis  the  block  diagram  was  reduced  to  standard 
form.  Note  that  an  all  sampled  data  system  is  assumed. 


where  G  = 
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The  characteristic  equation  of  the  system  then  becomes 


1  -  G,H  +  G.G0  H  =  0 
J.  I 


The  numerator  of  this  expression  is  the  only  part  of  interest  and 
after  combination  becomes 


z  +  (Kg  +  Kh  -  2)z  +  (1  -  Kg)  =  0 


r  +  1 


Transformation  to  r  -  space,  z  =  — - j  gives 


Kh  r  +  2Kg  r  +  (4  -  2Kg  -  Kh)  =  0  . 


For  stability  all  poles  must  be  in  the  left  half-plane,  so  all  coefficients 
above  must  have  the  same  sign.  For  positive  K  the  first  two  coefficients 

are  always  positive,  so  from  the  r°  coefficient 


K  < 


4 

2g  +  h 


This  stability  requirement  is  plotted  on  page  14  . 

The  k  must  now  be  related  to  the  radar  system  parameters.  From 
the  simulation: 


■  <SLf)  (&$) 

where  a/E  is  the  ratio  of  difference  to  sum  channel  voltages. 

For  antenna  1 

sin  (P^  -  P:P3)  sin  (Pxp  +  P1P3) 

sin  (P2^  “  P2P3)  sin  (P2^  +  P2P3) 

sin  (Pxp  -  PXP3)  sin  (P^  +  PjPj) 

"  “in  (PjP  -  P2P3)+  sin  (P ‘p  +  P2P3) 


1) 

AO) 

So) 
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,»v  i9/;rr*r&r  y 


, 'T \*  ' '■v^viTT'  — t  - ; »■  v’ ^  -  •*  fyr..v 


■’wi'^r  *— *  ^VT77i^pf(^yK-^y^5K^- v*-^  ',  •  - 'rf >» o-^ rp3^r<?j 1  '^TJTVmT-^R -T) *z*~rzi *^ •  ^; ‘ x » t * *■?? sip 


«s»7,Ttrow 


where 


By 

reduces  to 


„  Nits 

pisT 

p  *sa 

2  \ 

p  =  -*_ 
r3  2Ns 

N  *  number  of  elements  per  row  or  column 
s  =  element  spacing  for  rectangular  array 
X  =  wavelength. 

the  use  of  trigonometric  identities  the  £  using  the  above 


fi}  - cot  (a)  tan  M 


2<e) 

and  for  the  parameters  used  in  this  case 


=  45.8  tan  (1.8  p)  «  82.5  p 

The  small  angle  approximation  may  be  justified  since  no  discontinuity 
occurs  in  the  region  of  interest,  and  a  good  straight  line  fit  is 
possible.  Thus 


?“  kl? 


kx  «  82.5 


and  for  proper  interpretation  of  the  measurement 


SLP  =  k2  « 


This  gives  the  general  and  specific  expressions  for  the  transfer 
function: 


K  °  * 


*1 


82.5 


,  k2  a  0.01208 
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2)  For  Antenna  2  a  closed  form  expression  could  not 
easily  be  derived,  so  the  curve  was  evaluated  numerically  and  the 
results  shown  in  Figure  B-l.  The  curve  for  Antenna  I  is  also  shown  for 
comparison.  An  interesting  and  very  practical  problem  arises  when  the 
value  for  1«2  is  selected.  If  the  slope  (k^)  is  measured  at  p  =  0,  then 

=  0.02360  radian/volt/volt.  However,  as  indicated  on  the  figure, 

this  may  not  be  as  good  as  the  value  sometimes  used  by  passing  a  line 
through  the  origin  and  the  1-degree  point  and  using  this  slope,  in  this 
case,  k2  =  0.01630  radian/volt/volt. 

Several  simulation  runs  were  made  and  the  stability  criteria 
presented  in  this  case  verified. 


